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TECHNICAL PUBLICATION 


THE GLOBAL LAND-OCEAN TEMPERATURE INDEX IN RELATION TO SUNSPOT 
NUMBER, THE ATLANTIC MULTIDECADAL OSCILLATION INDEX, THE MAUNA LOA 
ATMOSPHERIC CONCENTRATION OF C0 2 , AND ANTHROPOGENIC 

CARBON EMISSIONS 


1. INTRODUCTION 


Global warming/climate change has been a subject of scientific interest since the early 19th 
century. 1-7 In particular, increases in the atmospheric concentration of carbon dioxide (C0 2 ) have 
long been thought to account for Earth’s increased warming, 8-16 although the lack of a dependable 
set of observational data was apparent as late as the mid 1950s. 17-20 However, beginning in the late 
1950s, being associated with the International Geophysical Year, the opportunity arose to begin accu- 
rate continuous monitoring of the Earth’s atmospheric concentration of C0 2 . 21 ^ 26 Consequently, it 
is now well established that the atmospheric concentration of C0 2 , while varying seasonally within 
any particular year, has steadily increased over time. 27-29 Associated with this rising trend in the 
atmospheric concentration of CO? is a rising trend in the surface-air and sea-surface temperatures 
(SSTs). 30 - 36 

This Technical Publication (TP) examines the statistical relationships between 10-year mov- 
ing averages (10-yma) of the Global Land-Ocean Temperature Index (GLOTI), sunspot number 
(SSN), the Atlantic Multidecadal Oscillation (AMO) index, and the Mauna Loa C0 2 (MLC02) 
index for the common interval 1964—2006, where the 10-yma values are used to indicate trends in the 
data. Scatter plots using the 10-yma values between GLOTI and each of the other parameters are 
determined, both as single-variate and multivariate fits. Scatter plots are also determined for MLC02 
using single-variate and bivariate (BV) fits, based on the GLOTI alone and the GLOTI in combina- 
tion with the AMO index. On the basis of the inferred preferential fits for MLC02, estimates for 
MLC02 are determined for the interval 1885-1964, thereby yielding an estimate of the preindustrial 
level of atmospheric concentration of C0 2 . Lastly, 10-yma values of MLC02 are compared against 
10-yma estimates of the total carbon emissions (TCE) 37 to determine the likelihood that manmade 
sources of carbon emissions are indeed responsible for the recent warming now being experienced. 
(Parametric values used in this TP are those available prior to the end of 2012.) 


2. RESULTS 


Figure 1 displays the variation of the annual (thin line) and 10-yma (thick line) values of 
(a) GLOTI, (b) SSN, (c) AMO, and (d) MLC02 for the common interval 1959-2011. Individual 
sunspot cycles are identified in subpanel (b). Inspection of figure 1 suggests that during the common 
interval 1959-2011, the GLOTI, AMO, and MLC02 are all trending upwards, while SSN is trend- 
ing downwards. Furthermore, because the 10-yma value of MLC02 remains above its annual (same 
year) average during this common interval, it is apparent that the atmospheric concentration of C0 2 
as measured at Mauna Loa, Hawaii, is presently increasing at an accelerated rate, as placement of 
a straight-edge along the curve clearly shows. 

The GLOTI is a measure of the global land-ocean temperature relative to the base period 
of 1951-1980, where the data are taken from the Global Historical Climate Network, version 3, 
using elimination of outliers and homogeneity adjustment. The data, as used in figure 1(a), are the 
January-December averages, available online at <http://data.giss.nasa.gov/gistemp>. 33 ’ 38 

The SSN provides a measure of the strength of solar activity associated with the vari- 
ation of solar irradiance over the solar cycle. 39 ^ 2 Annual values of SSN are available online at 
<ftp://ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SUNSPOT_NUMBERS/INTERNATIONAL>. 

The AMO is a temperature oscillation in the detrended SST of the North Atlantic Ocean 
(0-70° N.) having a period of about 65-70 years that fluctuates between warm (positive) and cool 
(negative) phases. 43 ’ 44 The SST pattern has been suggested to be linked to variations in the strength 
of the Atlantic thermohaline circulation (THC), a density-driven, global circulation pattern that 
involves the movement of warm equatorial surface waters to higher latitudes and the subsequent 
cooling and sinking of these waters in the deep ocean 45 In particular, the warm phase of the AMO 
represents intervals of faster THC, while the cold phase represents intervals of slower THC. From 
figure 1(c), one surmises that the AMO was reflective of the cool phase between about 1964 and 
1994, but now is in the midst of a warm phase (since 1995) that began its rise to more positive values 
about 1975 (based on 10-yma values). The warm phase is expected to continue for at least another 
decade or more. Monthly values of the AMO index are available online at <http://www.esrl.noaa.gov/ 
psd/data/correlation/amon.us.long.data>. 

The MLC02 is a measure of the atmospheric concentration of C0 2 as measured at the 
Mauna Loa Observatory on the Big Island of Hawaii, located on the northern slope of the volcano 
Mauna Loa at an elevation of 3,400 m above sea level and 800 m below its summit at 19.5° N. and 
155.6° W. 24 ’ 46 Annual averages of the MLC02 measurements, expressed in ppm, are available online 
for the interval 1959-2011 at <ftp://ftp.cmdl.noaa.gov/ccg/co2/trends/co2_annmean_mlo.txt>. 
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Figure 1. Variation of annual (thin line) and 10-yma (thick line) values of (a) GLOTI, 
(b) SSN, (c) AMO, and (d) MLC02 for the common interval 1959-2011. 
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Figure 2 depicts the scatter plots of the 10-yma values of GLOTI versus (a) SSN, (b) AMO, 
and (c) MLC02 for the interval 1964—2006. The diagonal lines in the subpanels represent the inferred 
preferential linear regressions. The small arrows in each subpanel (and in succeeding figures) indicate 
the last available 10-yma values for the parameters (i.e. , for the year 2006). The 10-yma values in each 
subpanel generally track smoothly from lower to higher GLOTI values over the interval 1964—2006. 
Given in each subpanel are the inferred regression equation (y), the coefficient of correlation (r), 
the coefficient of determination (r 2 ) (a measure of the amount of variance explained by the inferred 
regression), the standard error of estimate ( se ), and the confidence level (cl) of the inferred regres- 
sion, where cl >90% indicates a marginally statistically important result, cl> 95% indicates a statisti- 
cally important result, and cl >99% indicates a highly statistically important result. 




(b) AMO (10-yma) 



Figure 2. Scatter plots of the 10-yma values of GLOTI versus (a) SSN, 
(b) AMO, and (c) MLC02 for the interval 1964-2006. 
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For GLOTI versus SSN, the inferred regression based on the limited time interval 1964—2006 
is one found to be only of marginal statistical importance, having cl> 90%, se = 0.2 °C, and r 2 = 0.072, 
meaning that the inferred regression explains only about 7% of the variance in the 10-yma values of 
GLOTI during the interval 1964—2006. Also, the inferred regression is found to be one that varies 
inversely (r = - 0.27) rather than directly, in contrast to that found when using the much longer inter- 
val 1885-2006, one that clearly shows an upward-trending SSN. 47 

For GLOTI versus AMO, the inferred regression (r = 0.87) is one that is highly statistically 
important, having cl >99. 9%, sc = 0.09 °C, and r 2 = 0.762, meaning that the inferred regression 
explains about 76% of the variance in the 10-yma values of GLOTI. Current values of AMO and 
GLOTI are positive, indicative of the warm phase and the occurrence of a warm anomaly, respec- 
tively. The implication then is that the continuing warm phase as indicated by the AMO suggests 
a continuing positive (warm) global land-ocean temperature anomaly for the near term foreseeable 
future. 


For GLOTI versus MLC02, the inferred regression (r = 0.99) is the strongest of all the 
inferred regressions, having cl> 99.9%, .sc = 0.06 °C, and r 2 = 0.988, meaning that the inferred regres- 
sion explains nearly 99% of the variance in the 10-yma values of GLOTI for the interval 1964—2006. 
The implication then is that the continuing unabated increase in the atmospheric concentration of 
C0 2 as measured at Mauna Loa strongly suggests a continuing increase in the positive (warm) global 
land-ocean temperature for the long term foreseeable future. 
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Figure 3 shows the scatter plots of the 10-yma values of GLOTI versus selected BV fits for 
GLOTI, including (a) BV1, based on SSN and AMO, (b) BV2, based on SSN and MLC02, and 
(c) BV3, based on AMO and MLC02. The BV fit BV1, having R yU = 0.895 and S yU = 0.086 °C, 
represents a slight improvement over using AMO alone, while the BV fit BV2, having i? 12 = 0.994 
and 5’ v12 = 0.021 °C, and especially BV fit BV3, having i? vl2 = 0.999 and 5' vl2 = 0.01 °C, represent sig- 
nificant improvements over using MLC02 alone (by virtue of their inferred reduced standard error 
of estimates). (In figure 3, subscripts 1 and 2 refer to parameters 1 and 2, where 1 is SSN and 2 is 
AMO in BV1; SSN and MLC02 in BV2; and AMO and MLC02 in BV3.) 
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Figure 3. Scatter plots of the 10-yma values of GLOTI versus selected BV fits for GLOTI, 
including (a) BV1, based on SSN and AMO, (b) BV2, based on SSN and 
MLC02, and (c) BV3, based on AMO and MLC02. 


Figure 4 displays the scatter plot of the 10-yma values of GLOTI versus a trivariate (TV) fit 
for GLOTI (TV1), based on SSN, AMO, and MLC02. Essentially, the TV1 fit is insignificantly dif- 
ferent from that of the BV3 fit (i.e. , the r and se are essentially the same for both fits). 
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Figure 4. Scatter plot of the 10-yma values of GLOTI versus a 
TV fit for GLOTI (TV1), based on SSN, AMO, and 
MLC02. 
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Figure 5 depicts the scatter plots of the 10-yma values of MLC02 versus (a) GLOTI and 
(b) BV4, the BV ht based on both GLOTI and AMO. Of the two fits, the one between MLC02 
and GLOTI alone is the stronger, having r = 0.994 and se = 2 ppm, as compared to i ? vl 2 = 0.992 and 
5^22 = 2.4 ppm for BV4. The diagonal line in figure 5(a) is the inferred preferential linear regression 
line for MLC02 versus GLOTI. The regression equations and other statistical parameters are identi- 
fied in the subpanels. 
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(b) MLC02 (GLOTI, AMO)=BV4 


Figure 5. Scatter plots of the 10-yma values of MLC02 versus (a) GLOTI and (b) BV4, 
the BV fit based on GLOTI and AMO. 
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Figure 6 shows the estimated and observed 10-yma values of MLC02 for the interval 1885— 
2006 based on the inferred regressions shown in figure 5. Thus, presuming the validity of the inferred 
regressions, one finds that, while both fits essentially mimic the observed 10-yma values of MLC02 
(the filled circles) for the interval 1970-2006, slight differences are apparent in the values prior to 
1970. In particular, BV4 (the thick line) suggests slightly lower estimated values (by about 10 ppm) 
for the 10-yma values of MLC02 during the interval 1885-1964, as compared to using the linear 
fit based on GLOTI alone. Both extrapolations suggest temporary slight increases in atmospheric 
concentration of C0 2 near 1900 and the early 1940s, each increase lasting about 20 years in length. 
Values since about 1975 have been the highest on record and continue to increase with the passage of 
time. 



Year 


Figure 6. Estimated and observed 10-yma values of MLC02 for the interval 1885-2006, 
based on the inferred regressions shown in figure 5. 
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Figure 7 compares the estimated (based on BV4) and observed 10-yma values of MLC02 
for the interval 1885-2006 with exponentials based on year (t). Fit jq is based on the exponential 
determined using the 10-yma values from BV4 for t- 1890 and /= 1964, and ht y 2 is based on the 
exponential determined using 10-yma values from BV4 for t- 1964 and t = 2006. For 1890, 1964, 
and 2006, the 10-yma values from BV4 equal 285.25, 320.08, and 380.39 ppm, respectively. The hts 
jq and y 2 are provided merely to demonstrate the approximate fitting to the observed and estimated 
MLC02 based on using the BV4 fit and the apparent change that seems to take place in the trend- 
ing of MLC02 values relative to the exponential fits sometime around 1970, when the observed and 
estimated MLC02 values appear to have transitioned from one following more closely the y l fit to 
one following more closely the y 2 fit. (Actual 10-yma values of MLC02 are now found to be above 
the estimated values given by the y 2 exponential fit, indicating that higher values of MLC02 should 
be expected sooner rather than later, as given by the y 2 exponential fit.) 
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Figure 7. Comparison of the estimated (based on BV4) and observed 10-yma values 
of MLC02 for the interval 1885-2006 with exponentials based on year t. 


The expected 10-yma value of MLC02 about 1880 is of interest because it is representative of 
the preindustrial level of the atmospheric concentration of C0 2 . Now, the record of GLOTI yearly 
values begins in 1880, having a value of -0.32 °C. Consequently, the first available 10-yma value 
of GLOTI occurs in 1885. However, the record of AMO yearly values begins in 1856, so that the 
10-yma value of AMO is known for 1880 (=0.075 °C). For the interval 1885-2006, the residual (i.e., 
the difference between the yearly value and the corresponding 10-yma value) for GLOTI is about 
±0.1 °C, a range that captures about 80% of the residuals. 47 Hence, one can estimate the 10-yma 
value for GLOTI for the year 1880 to be about -0.32 ±0.1 °C. Using this value in conjunction with 
the 10-yma AMO value for 1880 in BV4 allows one to estimate the 10-yma value of MLC02 for 
1880, a value computed to be about 283.6 ±2.4 ppm, very close to that determined from polar ice 
cores for the preindustrial level (i.e., about 280 ± 5 ppmv; see Neftel et al. 48 ). 

As previously noted, the atmospheric concentration of C0 2 as measured at Mauna 
Loa appears to be not only increasing over time, but also to be increasing at an accelerated rate. 
Figure 8 plots the (a) first and (b) second differences of the 10-yma values of MLC02 for the interval 
1964—2006. From figure 8(a), one finds that the average 10-yma value of MLC02 is growing at the 
rate of about 1 .47 ±0.31 ppm per year (i.e., the mean ± 1 standard deviation (^interval). However, 
as can be seen from figure 8(a), the average rate inadequately describes the actual year-to-year first 
difference values. The rate of growth in the 10-yma value of MLC02 was always below 1.5 ppm 
prior to 1980 and, except for the brief intervals 1987-1989 and 1991-1992, it always has been equal 
to 1.5 ppm or greater after 1980, measuring 2.03 ppm in 2005 (the last available first difference value 
and the largest rate of growth ever measured). Overall, the first difference values are better described 
using the linear fit given by y = -46.207 + 0.024/, where y is the expected first difference and t is the 
year, rather than using the mean and sd. The linear fit is found to explain about 88% of the variance 
in the first difference values. For 2006, one expects the first difference to be about 1 .99 ±0.3 ppm (the 
±1 se prediction interval), thereby inferring that the expected 10-yma of MLC02 for 2007 should be 
about 381.6 + 1.99 ±0.3 ppm or about 383.59 ±0.3 ppm. Presuming the validity of the fit, one esti- 
mates the 10-yma value of MLC02 for 2015 and 2026 to be >400 ppm and >425 ppm, respectively. If 
these values prove true, then, from figure 2(c), one expects the 10-yma value of GLOTI to measure 
about 0.72 ±0.06 °C in 2015 and about 0.97 ±0.06 °C in 2026. (It should be noted that the 10-yma of 
MLC02 for 2007 is now known to measure 383.67 ppm, thus confirming the prediction above that 
it would measure about 383.59 ±0.3 ppm.) 
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Figure 8. The (a) first and (b) second differences of the 10-yma values of MLC02 
for the interval 1964-2006. 


Figure 8(b) displays the variation of the second difference of the 10-yma values of MLC02 
(i.e., the first difference of the first difference). It has a mean of 0.027 ±0.047 ppm per year. Runs- 
testing 49 of the second difference values strongly suggests that they are distributed randomly, hav- 
ing a normal deviate for the sample (z) = -0.467, based on the observed number of values above the 
median (na = 21), the number of values below the median (nb- 20), and the number of runs above 
the median (nra= 10), where the median second difference value is 0.03. 


12 


Recently, Boden et al. reported on the yearly variation of global fossil fuel C0 2 emissions. 37 
They noted that approximately 356 billion metric tons of carbon have been released into the atmo- 
sphere through the combustion of fossil fuels and cement production since 1751, with half of the 
fossil fuel C0 2 emissions having occurred since the mid 1980s. Furthermore, they noted that liquid 
and solid fuels accounted for about 76.6% of the emissions in 2009, while gas fuels (e.g., natural gas) 
and cement production accounted for 17.9% and 4.7% of the emissions, respectively. 

Figure 9 plots the annual and 10-yma values of the TCE estimates for the interval 1 959— 
2009. For the overall interval 1964-2004, the trend in TCE has grown at the average rate of about 
118 million metric tons per year, although the intervals prior to 1975 and after about 1997 have 
had considerably higher average rates of growth, measuring about 157 and 169 million metric tons 
per year, respectively. Extrapolation of the 10-yma trend line beyond 2004 using the higher rate of 
growth in TCE suggests that TCE will measure about 10,000 million metric tons around the year 2017 
± 2 years, and it will measure about 12,000 million metric tons around 2029 ± 2 years. (The rightmost 
vertical axis gives the total CCE emission in units of million metric tons.) 
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Figure 9. Annual and 10-yma values of the TCE estimates for the interval 1959-2009 
(taken from Boden, Marland, and Andres, 20 12). 37 
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Figure 10 plots the 10-yma values of MLC02 versus TCE, assuming a zero time lag. Clearly, 
as TCE increases, MLC02 increases, strongly suggesting that the two parameters are highly statisti- 
cally correlated, which they appear to be (r = 0.985 and cl» 99.9%). The lower-leftmost tilled circle 
corresponds to the year 1964, and the upper-rightmost tilled circle corresponds to the year 2004, with 
the intervening filled circles corresponding to the intervening sequential years. From the inferred 
regression (y ), one finds that a 10-yma value of TCE= 10,000 million metric tons corresponds to 
a 10-yma value of about MLC02 = 408.75 ± 3.79 ppm (the ±1 se prediction interval), while a 10-yma 
value of TCE= 12,000 million metric tons corresponds to about MLC02 = 436.31 ± 3.79 ppm. Simi- 
larly, from the inferred regression (x), one finds that a 10-yma value of MLC02 of 400 ppm cor- 
responds to a 10-yma value of about TCE = 9,246 ±21 8 million metric tons, while 10-yma values of 
MLC02 = 425 ppm and 440 ppm correspond to about TCE= 11,006 ±218 and 12,062 ±21 8 million 
metric tons, respectively. 



Figure 10. Scatter plot of the 10-yma values of MLC02 versus TCE and TCE versus MLC02. 


As previously noted from figure 9, extrapolation of the 10-yma trend line beyond 2004 using 
the higher rate of growth in TCE suggests that TCE will measure about 10,000 million metric tons 
around the year 2017 ±2 years, and it will measure about 12,000 million metric tons around 2029 
±2 years, presuming the validity of the extrapolation. Furthermore, from figure 2(c), given the 
inferred 10-yma values of MLC02 corresponding to 10-yma values of TCE equal to 10,000 metric 
tons and 12,000 million metric tons, being equal to about 409 ppm and 440 ppm, respectively, one 
estimates 10-yma values of GLOTI equal to about 0.82 ±0.06 °C around 2017 ±2 years and about 
1.13 ±0.06 °C around 2029 ±2 years, respectively. A faster (slower) rate of growth in the 10-yma 
values of TCE implies higher (lower) 10-yma values of MLC02 and, consequently, higher (lower) 
10-yma values GLOTI. 
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3. DISCUSSION AND CONCLUSION 


This TP of the common interval 1964-2006 has revealed a close statistical coupling between 
the global land-ocean temperature anomaly as described by GLOTI and both the phase of the AMO 
and the increasing atmospheric concentration of C0 2 as measured at Mauna Loa. On the basis of 
10-yma values, the GLOTI is found to be well described using the BV ht BV3 = -2.77 + 0.224552AMO 
+ 0.008615MLC02, having 7^^ = 0.9986 and 5 , vl2 = 0.01 °C, where BV3 represents the 10-yma of 
value of the GLOTI, and the subscripts 1 and 2 refer to the 1 0-yma parametric values of the AMO and 
MLC02 indices, respectively. Likewise, the BV ht BV4 (=322.309 + 1 15. 027GLOTI-25.038AMO) 
is found to provide a very close ht for the observed 10-yma values of MLC02, having i? rl2 = 0.992 
and S yl2~ 2 .42 ppm. 

From an extrapolation of the BV4 ht backwards in time, prior to 1964, one hnds evidence for 
two temporary slight increases in the atmospheric concentration of C0 2 , each of about 20 years dura- 
tion — the hrst about 1900, having a local peak of about 296 ±2.4 ppm, and the second in the early 
1940s, having a local peak of about 322 ± 2.4 ppm. While the interval prior to about 1970 appears to 
be better approximated using the exponential log y j = 1 .856 + 0.002/, where y l is the 10-yma value of 
MLC02 and / is the year, the interval after 1970 appears to be better described using the exponential 
y 2 =-2.085 + 0.004/. However, current 10-yma values of MLC02 now appear to be rising faster than 
described by the post-1970 exponential. Using the BV4 ht, one estimates the preindustrial value for 
the atmospheric concentration of C0 2 to be about 283.6 ±2.4 ppm, very close to that determined 
from polar ice cores (i.e., about 280 ± 5 ppmv). 

The 10-yma value of MLC02 measured about 382 ppm in 2006. On average, the 10-yma 
value of MLC02 is found to be increasing at the approximate rate of about 1.47 ±0.31 ppm per 
year. While true, the average rate of growth in the 10-yma values of MLC02 does not appear to 
adequately describe the actual year-to-year increases in the 10-yma values. Instead, the rate of yearly 
increase is better described using the linear ht y = -46.207 + 0.024/, where y is the hrst difference in 
consecutive year-to-year 10-yma values of MLC02 and / is the year. Hence, the 10-yma values of 
MLC02 appear to be increasing at an accelerated rate, such that the 10-yma value of MLC02 is 
expected to measure >400 ppm about the middle-to-latter half of the present decade and >425 ppm 
about the middle-to-latter half of the following decade. Likewise, because of the inferred close asso- 
ciation between the 10-yma values of GLOTI and MLC02, one expects 10-yma values of GLOTI to 
continue to increase over time, measuring about 0.72 ± 0.06 °C and 0.97 ± 0.06 °C, respectively, when 
10-yma values of MLC02 measure >400 ppm and >425 ppm. 

The rise in 10-yma values of MLC02 appears to be directly linked to the increase in the 
10-yma values of the TCE as deduced from the yearly values of global fossil fuel C0 2 emissions 
determined by Boden et al. 37 In particular, the 10-yma values of MLC02 are found to increase in 
step with increases in the 10-yma values of TCE, with the inferred correlation having r = 0.985 and 
cl» 99.9%. From the inferred regression between the 10-yma values of MLC02 and TCE, one finds 
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that a 10-yma value of TCE equal to about 10,000 million metric tons suggests a 10-yma value of 
MLC02 >400 ppm and GLOTI >0.76 °C, while a 10-yma value of TCE equal to about 12,000 mil- 
lion metric tons suggests a 10-yma value of MLC02 >440 ppm and GLOTI >1.07 °C. 
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